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1 Introduction 

This paper deals with phenomena of phase transitions, of first and second order, in binary 
systems (cf., e.g., \2U \ 123"]). In a first order phase transition phenomenon, as in the solid- 
liquid or liquid-vapor phase change, the phase transition occurs at a critical temperature, say 
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$ c : if the absolute temperature $ in the body is strictly greater than the critical temperature i? c , 
then the minimum of the energy potential is attained in one of the pure phases, while if $ < 1? C 
the minimum is attained in the other phase. In the case when i? = i? c the energy potential has 
two minima attained for the two phases, that is phase change may occur. On the other hand, in 
the case of second order phase transitions, the system behaves differently provided i? is greater 
or less than the critical temperature i? c . Indeed, for high temperatures the energy potential has 
only one minimum, while for # < # c two minima are attained with the same values. This second 
behaviour is characteristic, for instance, of some solid-solid phase transitions, ferromagnetism, 
and superconductivity. 

We are going to investigate a model describing these phenomena by use of phase- field theories, 
in terms of temperature i? and a phase parameter X, that includes the effects of micro-motions 
and micro- forces responsible for the phase transition (cf. [22] and [25] ). Indeed, it is known that 
phase transitions are caused by changes occurring at a microscopic level in the (atomic and/or 
crystal) structure of the system. These changes are the effects of micro- forces and motions, which 
have to be included in the balance of the energy of the whole system, even if we are providing 
a macroscopic description of the phenomenon. We follow the suggestion of M. Fremond, who 
proposed a new balance law on phase transitions (cf. [22] and [16] ) . Let us quote a fairly similar 
theory devised by Gurtin [25] for Ginzburg-Landau and Cahn-Hilliard equations. 

Hence, we combine this theory with a model, recently introduced, based on a reduced energy 
balance equation, subsequently termed as entropy equation since entropy is involved. We mainly 
refer to [8] and especially [5| for the derivation of the model and related analytical results. 
Note that in [5] also thermal memory effects are accounted for (according to the theory of 
|26j). while in the present contribution we are neglecting them. The main advantage of the 
model itself is that, once the problem is solved in a suitable sense, one can obtain directly 
the positivity of the temperature, mainly due to the presence of a logarithmic nonlinearity in 
the resulting system of partial differential equations. This avoids the application of maximum 
principle arguments, which are difficult to set in a number of interesting situations. This fact 
is pointed out in [5], where the model has been introduced in a general setting and a global 
existence result is proved for a weak formulation. For the sake of completeness, we quote a 
recent contribution [9] , where a more general (non-smooth) relation between the entropy and the 
absolute temperature is considered. Also in this case the physical constraint on the positivity of 
the absolute temperature in ensured by the model itself. In [9] the model is recovered by writing 
the first principle in a dual formulation (in the sense of the convex analysis) using the entropy in 
place of the temperature as state variable. The authors of [9j prove existence of solutions (t?, X) 
and provide a characterization of the long-time behaviour of the solution trajectories (see also 
[U[TU]). However, uniqueness is still an open problem (both in [5] and in [S]): this is mainly due 
to the lack of regularity on the •& component of the solution. In [6] [7] a suitable choice (far from 
the present approach) of the heat flux law leads to a linear operator acting on the temperature, 
which is of some help in showing existence, uniqueness, and regularity of the solution. Also the 
large time behaviour and the w-limit set are investigated in [7J. 

The model. In the actual contribution, we mainly proceed according to [5]. Indeed, we con- 
sider a two-phases system, located in a smooth and bounded domain Q C M 3 , and look at its 
evolution during a finite time interval (0, T). We denote by T the boundary dQ. The thermo- 
mechanical equilibrium of the system is described in terms of state variables and governed by 
the free energy, while the dynamics ensues from the presence of a pseudo-potential of dissipation 
(depending on dissipative variables). Let us point out that the properties of the pseudo-potential 
of dissipation ensure thermodynamical consistency of the model. 

We do not consider mechanical effects, so that the variables of the system are just the 
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absolute temperature $ and a phase parameter X, related to the proportion of one phase with 
respect to the other. In general, X attains its physical admissible values in a range [X*,X*] (e.g., 
X € [0, 1]) and this physical constraint has to be ensured by the model itself. Hence, we derive 
the equations of the system by thermomechanical laws. More precisely, we use the approach 
followed by Fremond (yielding the evolution of the phase parameter), and the first and second 
principle of Thermodynamics (from which we recover a balance equation ruling the evolution of 
the temperature). 

We first discuss, with some details, the derivation of the evolution equation for the phase 
parameter X. It is known that the phase transition can be described as a change in the order 
structure of the thermomechanical system we are considering, so that the phase parameter X 
can be interpreted also as an "order parameter". More precisely, below a critical temperature 
it is observed that the structure order of many materials is greater than above. For instance, 
in the solid-liquid phase transition the solid phase has a greater order structure due to the 
crystal symmetry group. Analogously, we may think to ferromagnetism where the magnetic 
moments are aligned below the Curie temperature. Obviously, order change in the structure of 
the system occurs at a microscopic level. However, the parameter X is a macroscopic parameter 
whose evolution is governed by a balance law on the order structure, responsible for the phase 
transition at a microscopic level. Thus, the evolution of the order parameter X can be derived 
from the balance conditions 

£-divH = infix(0,T), H n = onrx(0,T) (1.1) 

n being the outward normal to the boundary T, where the scalar quantity B can be interpreted 
as an internal order structure density per units of concentration X, and H represents an order 
structure flux vector. Note that here we are assuming that no external action is provided to the 
system. 

As far as the description of the evolution of the temperature is concerned, we use a balance 
law in which higher order dissipative contributions are neglected by means of the small pertur- 
bations assumption (cf. [23]). The resulting equation rules the evolution of the entropy s of the 
system in terms of the entropy flux Q and an external source R, possibly depending on the state 
field, that is 

s t + divQ = i? inOx(0,T). (1.2) 

Then, (II. 2D is combined with boundary conditions on the entropy flux, e.g., if no flux is assumed 
through the boundary, then one states that Q • n = on T x (0, T). Another possibility, actually 
followed by us in the analysis, is to prescribe the value of the temperature on the boundary. 

Referring to [5] and [9], we prefer not to detail here the derivation of the model. However, 
for the sake of completeness, we just recall that (jl.2p can be obtained, in the framework of small 
perturbations assumptions, dividing by •& > the energy balance 

e t + div q = r + BX t + H-VX t . (1.3) 

In (II. 3D e denotes the internal energy of the system, q = $Q represents the heat flux, r stands 
for an external source, while BXt + H • VXt is internal and comes from the order structure. We 
will add some comments below on the relation between the two equations (|1.2|) and (|1.3p . 

Free energy. We specify the involved physical quantities with the help of two energy func- 
tionals: the free energy vf, depending on the state variables and accounting for the thermo- 
mechanical equilibrium of the system, and the pseudo-potential of dissipation 5> (see [28|). 
defined for dissipative variables and responsible of the evolution of the system. More precisely, 
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we consider as state variables the absolute temperature the order parameter X, and its gra- 
dient VX. It is known from Thermodynamics that the free energy is a concave function with 
respect to the temperature, while there are not constraints concerning the dependence on the 
other variables. In the present contribution, we choose the functional *I> of the following form 

¥(0,X,VX) := -y?? 2 + F{X)$ C + G{X)$ + ^|VX| 2 (1.4) 

the constants cq, v being positive and i? c > representing the critical value of the temperature 
for the phase transition. Note that the purely caloric part in the free energy — (co/2)i9 2 is, in 
fact, concave with respect to the temperature. Then, the functions F and G characterize the 
behaviour of the phase transition. For instance, in a first order phase transition, as vapor-liquid 
or liquid-solid, they can be prescribed as follows 

. A A „ , . A Zi/X. A . . 

F(X) = , G(X) = + — 1.5) 

K J 4 3 > \ J 4 3 t- 2 \ J 

while in a second order phase transition, as for superconductivity or ferromagnetism, F and G 
can be written as 

F(X) = -- Y , G(X) = -. (1.6) 

Let us remark from the beginning that both the cases (jl.5p and (jl.6p comply with our general 
assumptions (j2.3p ~ fl2.4p provided you take X* < and X* > 1. 

Admissible values for the phase variable. The physical constraint on X, that is X* < X < 
X*, is not a priori guaranteed by the choice of the free energy (jl.4p : in particular, the functions F 
and G prescribed in (jl.5p and fll.6j) are smooth on the whole of R (as in [3|). A different possible 
choice for ^ introduced in the literature (cf., e.g., [221 15]) and accounting for this constraint is 

*(tf,X,VX) = -|tf 2 + / [XtX] (X) +3"(X) + tfS(X) + ^|VX| 2 

where 5F, 9 are sufficiently smooth functions characterizing the phase transition. In this case 
the free energy is defined for any value of X but it is +oo if X [X*, X*] (while 2-/^ %*](X) = 
if X S [X*,X*]). In our approach, instead, we will show that the constraint on X is ensured by 
the evolution of the system, i.e., it will be proved that the equations of the system are somehow 
consistent as they yield X* < X < X*. The proof of this property of our model will be detailed 
in the sequel and relies on a maximum principle argument. 

Pseudo-potential of dissipation. Secondly, we introduce the pseudo-potential of dissipation 
$ (see |28j ) that depends on the dissipative phase variables Xt and Vi9. Let us just comment 
on the choice of these dissipative variables: Xt is related to microscopic transformations which 
are responsible for the phase transition, i.e. for the evolution of the order structure of the 
system, while V$ is concerned with the heat flux. For the sake of completeness, let us recall 
that the pseudo-potential of dissipation is non-negative, convex with respect to the dissipative 
variables, and it attains its minimum for a null dissipation, that is when (Xt, Vi9) = (0, 0). We 
prescribe 

d.(X i ,V^) = ^|X t | 2 + A|v^ (i.7) 
with and A denoting positive coefficients. 

Constitutive relations. Hence, constitutive relations can be written for B, H, s, Q. They 
are recovered from the free energy (for non dissipative contributions) and the pseudo-potential 
of dissipation (for dissipative parts). We have 

B= dx + dx t = ^ c (x) + (x) + (L8) 
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and 

H = — . = uVX (1.9) 

as well as 

s = ~ = c tf - G(X) (1.10) 

and 

Let us point out that the choice of the free energy (jl.4[) leads to a linear contribution for the 
temperature in (jl.lOp . This will yield sufficient regularity on the solution, from which we will 
be able to prove uniqueness. We also point out that the term — (co/2)$ 2 could be seen as a first 
order approximation of the following, well-known, form of the energy potential 

¥(0, ...) = _ Co tfiogtf+ ... 

used, e.g., in [5]. In this case, the entropy s would be related to the temperature i? through a 
logarithmic nonlinearity. Note that, on the contrary, a logarithmic nonlinearity forcing $ to be 
strictly positive is present in our expression (jl.lip for Q. 

From energy balance to our equation. The reader may be curious about the derivation of 
(|1.2|) from (II. 3p . Let us recall constitutive relations (|1.8j) - (|l.lip and the well known Helmoltz 
relation 

e = ^ + tfs. (1.12) 
By applying the chain rule in (|1.3f) , some terms cancel and one can rewrite the energy balance as 

fl („ + divQ- B )=(|i,^).(X„V*) (1.13) 



letting Q = q/i? and R = r/fi. We point out that the fact that is convex, l.s.c, proper, 
non-negative and it attains its minimum when (X^,V0) = (0,0) ensure that the right hand 
side of (jl . 131) is non-negative. As > 0, the Clausius-Duhem inequality 

s t + divQ- R>0. (1.14) 

complies with ()1.13j) . Hence, dividing (| 1 . 13[) by i? and neglecting the resulting contribution on 
the right hand side, equation (|1.2p follows. 



Nonlinearity in the source term. We stress that the entropy source R in (II. 2p is related to 
the source r appearing in (jl.3p by R = r/0. Thus, it seems reasonable to include in our analysis 
the possibility for R to depend on the temperature (and possibly to present some singularities). 
This is one of the features of our paper: we actually deal with entropy sources (positive or 
negative, as sinks) depending also on the temperature <&. Indeed, possible choices satisfying our 
assumptions are (cf. the later Remark 12 .2p 



R(x,t^) = ^^-R 2 (x,t) 
which would correspond to r(x,i,$) = (Ri(x,t)/$) — R2(x,t)$, or 



R(x, t, 0) = R 3 (x, t)tf - R 4 {x, t) 
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which can be viewed as a linearization of R around some equilibrium value of *&. In such cases, 
the possible data R\, R2 or R%, R4 are smooth enough and at least R\ should be non-negative 
throughout U x (0,T). 

System of PDEs and initial-boundary value problem. Now, combining constitutive 
relations (fTTH ]) - ([TTTT]l with (fTTTj) and (fL2|) leads to the following PDE system 

c $ t - G'(X)X t - AA log ■& = R{x, t, 0) (1.15) 
fiX t - uAX + F'(X)$ C + G'{X)ti = (1.16) 

which is addressed in Q := O x (0,7"). Then, concerning boundary conditions, prescribed 
in r x (0, T), we fix a Dirichlet condition for the temperature and a Neumann homogeneous 
condition for the phase parameter 

log?? = logtfr, d n X = 0. (1.17) 

Finally, initial conditions are set in Q 

0(0) = O , X(0) = X . (1.18) 

For the sake of simplicity, in the mathematical analysis performed in subsequent sections we will 
take the physical constants Co, A, fx, u, C all equal to 1. 

Let us point out in our model is that we are dealing with a non-smooth entropy source 
R(&) in (jl . 151) that is assumed to be increasing with respect to G (0, +00) up to Lipschitz 
perturbations. This fact turns out to be interesting both for analytical and modelling aspects. 
Indeed, we are able to treat a diffusive equation for the temperature with nonlinear and singular 
diffusion along with a non-smooth contribution in source term. Then, for modelling aspects, 
we can figure the entropy source R directly as r(0)/0 (compare (j 1 . 1 3 1) with (|1.3|) ). The second 
aspect we notice is that the physical property X £ [X*, X*] comes as a consequence of our analysis 
and then it is ensured by the evolution of the system itself and not required a priori. 

Under suitable assumptions on the data and on the regularity of the involved nonlinearities, 
we can prove existence of a solution in any time interval (0, T) for a weak formulation of our 
system (|1.15|) - (|1.18|) . Moreover, in a wide setting we are able to show that also the solution 
component is bounded by adapting a Moser type argument. Finally, if R(x, t, 0) is uniformly 
Lipschitz continuous with respect to 0, we prove uniqueness of the solution. 

Phase field model with special heat flux law. What is also interesting in our contribu- 
tion is that system (|1.15p - (|1.16p may be interpreted as a nonlinear Caginalp phase-field model 
(see [13]) with special heat flux law. Indeed, take the following energy functional 

1 wd 
T(0, X, VX) = -— + 0F(X) + G(X) + -^-|VX| 2 . (1.19) 

Hence, arguing as in the derivation of the Penrose-Fife model (cf. |29 } 112 } I2T]). we recover the 
entropy s, the internal energy e, and the phase field equation (for X) as 

<9T 1 „. , v .„ l9 

s = -m=-W- F{X) -2 m < L20) 

e = T + tfs = -i+G(X) (1.21) 
,X 1 + f = 0, ff = /J T (1.22) 
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where 5H /SX denotes the variational (Frechet) derivative of the functional H. Then, we can 
write the first principle neglecting microscopic motions and forces (here r represents a mere 
source term) 

ej + divq = r. (1.23) 

We now let 

q=-AVlogtf (1.24) 

and this special choice of the heat flux can be compared with the ones studied in the papers 
[141 115j . where Penrose-Fife models with special heat flux laws have been investigated. Then, 
combine (fL22j) and ([03]) with (TOI - tOTI) and (Ti~24"j) . so to obtain the following system 

u t - G'(X)X t - AA log u = -r (1.25) 
liX t -vAX + F'(X) + uG'(X) = (1.26) 

in which u has now the physical meaning of l/-d. Note that system fjl.25p — (|1.26j) is formally 
equivalent to (jl.l5j) - (j!.16j) (with cq = # c = 1 and — r in place of R). Thus, it turns out that 
our analysis applies to the Caginalp system (|1.25fl — fjl.26j) in which non-smooth heat sources 
depending on the temperature are admitted. 

Plan of the paper. We conclude the Introduction by giving an outline of the paper. In 
Section [2] we set our problem in a Sobolev spaces framework and make precise assumptions on 
the data and the involved functions, also stating main results. In Section we introduce an 
approximating problem and prove the existence of solutions for it. Next, Section H] is devoted 
to derive some a priori estimates and to pass to the limit as the approximating parameter goes 
to 0. A maximum principle is also checked to show that the constraint on X is always satisfied. 
Boundedness of i9 and uniqueness proofs are the subjects of the last two Sections [5] and [61 
respectively. 



2 Main results 

In this section, we carefully describe the problem we are going to deal with and state our results. 
First of all, we introduce the notation regarding the domain in which the evolution is considered. 
In the sequel, Q is a bounded open set in R 3 , whose boundary T is assumed to be of class C 2 . 
Moreover, d n is the (say, outward) normal derivative on T. Given a finite final time T, we set 
for convenience 

Q t := Q x (0, t) for every t G (0, T] and Q ■= Q T . (2.1) 

Next, we describe the structure of our system. We are given four constants X*, X* € M 

such that 

< tf* < 1 < i?* and X* < X* (2.2) 
and four functions F, G, (5, and it 

F, G : R -f R, 0:Qx(O, +oo) -»R, and vr : Q x R -► R 
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satisfying 

F,Gg C 2 (M), F is bounded from below and G is nonnegative (2-3) 

F',G'<0 in(-oo,X*), and F', G' > in(X*,+oo) (2.4) 

/3 is Lipschitz continuous in Q x [5, 1/5] for every 5 G (0, 1) (2-5) 

P,x , P,t, /3', and 7r are Caratheodory functions, with the notation (2-6) 

P, x (x,t,r) := V0(x,t,r), t (x, t, r) := d t (3(x, t, r), p'{x, t, r) := d r (3(x, t,r) (2.7) 
< P'{x,t,r) < /3i(r) 

for a.a. (x,t) E Q, every r£R, and some 0\ G C°(0,+oo) (2.8) 
\P, x (x,t,r)\ + \P,t(x,t,r)\ < Mp(l + \P(x,t,r)\) 

for a.a. (x,t) G Q, every r G (0, +oo), and some Mp G [0, +oo) (2-9) 

P(x, t,l) =0 for every (x, t) G Q. (2.10) 
|7r(x,t,r)| < A|r| +7T (x,t) 

for a.a. (x,t) G Q, every r£l, some A G [0, +oo), and some 7To G L 2 (Q). (2-H) 

Furthermore, we set for convenience 

P(x,t,r) := / P(x,t,s)ds for every (x,t,r) G Q X (0,+oo). (2.12) 



Then, observing that (by (j2.8j) ) P(x,t, •) is nondecreasing, it turns out that /? is nonnegative 
and convex with respect to the third variable. 

Remark 2.1. We note that the second inequality of (|2.8[) and (12.100 imply that 



(x,t,r)\ <P (r) :-- 



Pi(s) ds 



V (x,t,r) G Q x (0,+oo). (2.13) 



Therefore, we see by (|2.9p that even /? x and Pt satisfy an analogous inequality and infer that 
()2.5p follows from the other assumptions. We have written (12.50 in advance to give a meaning 
to the pointwise values of p. 

Remark 2.2. In the applications just the difference R = it — is interesting and we observe 
that the choice 

R (x t^i 

R(x,t,r) = U 2 ' ; +R 2 (x,t) withiii>0 (2.14) 

is included in our assumption. Indeed, we can take P(x,t,r) = R\(x,t) — R\(x,t) /r 2 for r > 
and n(x,t,r) = R±(x, t) +R%(x, t), so that R = it — and (|2.1(jp hold at the same time. However, 
we have to assume that R\ is Lipschitz continuous and that |Vi?i| + |<9t-Ri| < cR\ in Q for some 
constant c in order to fulfill (|2.9p . As Ri is bounded and nonnegative, (|2.8p holds as well. 
Moreover, we ask that R 2 G L 2 (Q) in view of ()2.1ip . A different choice for R is the following 

R(x,t,r) = R 3 (x,t)r - R 4 (x,t) (2.15) 

and can be obtained in our setting by simply taking = and tt = R. Here, we can assume 
R 3 G L°°(Q) and R A G L 2 (Q). 

Notation 2.3. Let / be a real interval and ^:Qx/— >Mbea Caratheodory function. We use 
the same symbol tp to denote the operator acting on measurable functions on Q as follows. If 
v : Q — > R is measurable 



tp(v) denotes the function (x, t) i— > ip(x, t, v(x, t)), (x,t)£Q. (2-16) 
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Note that ip{v) actually is measurable due to the Caratheodory assumption on tp. Similar 
definitions and symbols are used for functions depending on the space variable. Namely, if 
v : — > R is measurable 

tp(t,v) denotes the function x t— > ip(x, t, v(x)), (2-17) 

for a. a. t G (0, T). We obtain a time dependent operator. As a consequence, if v G L 2 (Q), the 
symbol ip(t,v(t)) denotes the measurable function x i— ► ip(x,t,v(x,t)), x G fi. Notation f|2. 16|) 
(|2.17p will be used, in particular, with (some of the functions listed below will be introduced 
later on) 

i> = ft ft ft, ft, fte, ft, ft, ft,,, ft,*, ft, vr. (2.18) 

Furthermore, we set 

H:=L 2 (Q), V'-H 1 ^), Vq := Hq(Q), W := {v £ H 2 (Q) : d n v = 0} . (2.19) 

We endow H, V, and W with their usual scalar products and norms and use a self-explaining 
notation, like || • ||y. For the sake of simplicity, the same symbol will be used both for a space 
and for any power of it. It is understood that H is embedded in V * in the usual way, i.e., 
(u, v) = (u, v) for every u G H and v G Vq, where ( • , ■ ) is the duality pairing between Vq and Vq 
and ( • , • ) is the inner product of H. 

Now, we list our assumptions on the boundary and initial data. We are given three functions 



fir, "&o, and Xq such that 

$ v £ L 2 {<d,T-H 1/2 (T))C\H l {Q,T-H~ l/2 {T)), 0* < r < 0* a.e. on T x (0, T) (2.20) 

00 G H, 0* < $ < 0* a.e. in O (2.21) 

X e V, X* < X < X* a.e. in (2.22) 

where 0*, 0*, X*, and X* are introduced in A2.2H . 

At this point, we are ready to state our problem. We look for a triplet (0, X,£) satisfying 
the regularity conditions and the equations listed below. 

G L°°(0, T; if), > a.e. in Q, and ln0 G £ 2 (0, T; V) (2.23) 

X G L 2 (0, T; W) (1 ii 1 (0, T; H) (2.24) 

G(X), F'(X), G'{X) G L 2 (Q) (2.25) 

£ G L 2 (Q) (2.26) 

^-G(X))6L 2 (0,T;^) (2.27) 

dt(d-G(X)) - Aln0 + £ = 7r(0) in L 2 (0,T;F *) and £ = ft0) (2.28) 

d t X- AX + F'(X) + G'(X)0 = a.e. in Q (2.29) 

In0 = hi0 r a.e. on T x (0, T) (2.30) 

(0 - G(X)) (0) = #o - G(X ) and X(0) = X . (2.31) 



Even though £ is a known function of 0, we refer to the triplet (0, X, £) instead of the pair (0, X) 
when we speak of a solution, just for convenience. Moreover, we note that (|2.23p and (|2.25p yield 
G'(X)0 G L 2 (0, T; L 1 ^)). However, we see by comparison in (|2~29jl that G"(X)0 G L 2 (Q). Next, 
we observe that G(Xq) makes sense in L°°(£l). Furthermore, we point out that the first condition 
in (|X3Tj) reduces to 0(0) = O whenever one knows that G(X) G C°([0, T]; F *). Actually, some 
additional smoothness for G(X) as well as for F'(X) and G'{X) surely holds if the nonlinearities 
satisfy some growth conditions, thanks to (|2.24|) . The same is trivially true whenever X is 
bounded, and our existence result stated below ensures such a property. Finally, (12. 24ft itself 
entails the homogeneous Neumann boundary condition for X (see (|2.19p ). Here is our first result. 
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Theorem 2.4. Let (|22 ]) - (|2TTT]) and (pi2"U]) - (f2~2"2"j) 6e fulfilled. Then, there exists a triplet 
satisfying the regularity requirements (j2.23[) ( f"2.27j) and solving problem (I2.28f ) — ( |2.31l) . 

Moreover, every solution X, £) fulfils the inequalities 

X*<X<X* a.e. in Q. (2.32) 

In particular, X is bounded. 

One can wonder whether the component $ of the solution is bounded as well. Actually, such 
a property holds whenever we reinforce the assumption on the structure of our system a little, 
namely 

7T G L q (Q) for some q > 5/2. (2.33) 
We can state the following result. 

Theorem 2.5. Assume (12.330 m addition to the hypotheses of Theorem \2.4\ Then, the compo- 
nent •& of any solution ($, to problem (j2.28j) (j2.31|) is bounded. 

Remark 2.6. From Theorem 12.41 it follows that G(X), F'(X), and G'(X) are smoother, namely 

G(X),F'(X), G'(X) G L°°(Q) n L 2 (0, T; V) n fl^O, T; £f) (2.34) 

due to (|2.3p and (|2.24p . In particular, we deduce that 

d t $ (£ L 2 (0,T;V *) (2.35) 

by comparison in (|2.27p . On the other hand, using the regularity conditions In t? G L 2 (0,T;V) 
and i? G L°°(Q) (see (l2~23|) and Theorem ESJ) , we see that W m = mi? m Vln^ for every m > 0, 
whence 

reI 2 (0,T ; y) for every m G (0,+co). (2.36) 

Moreover, as F'{X) and G'(X) are bounded, from (j!T24"]) and (f2T29|) we see that d t X - AX is 
bounded too. Hence, we have 

X e L p (0,T;W 2 ' p (n))nW^ p (0,T;L p (n)) for every p>l (2.37) 

thanks to the general theory of linear parabolic equations. 

Finally, we state our uniqueness result. As it often happens for doubly nonlinear problems, 
uniqueness for solutions cannot be proved unless more restrictive assumption on the structure 
are made. In particular, we cannot allow a singular (3 like (|2.14p . 

Theorem 2.7. Assume that the hypotheses of Theorem \2.J\ are fulfilled. Moreover, assume that 
R := 7r — (3 satisfies the uniform Lipschitz condition 

\R(x, t, r) — R(x, t,s)\ < \ft\r — s\ for every (x, t) G Q and r, s G (0, +oo) (2.38) 

for some \r > 0. Then, problem (|2.28| )-( f2.3ip has at most one solution ($, X, £) satisfying 
•& G L°°(Q). 

Corollary 2.8. If (I2.33H and (12.380 are fulfilled at the same time in addition to the hypotheses 
of Theorem \2.4\ then problem (|2.28|) - (|2.3ip has exactly one solution X,£) such that •& is 
bounded. Moreover, (]2.32|) . ()2.36p . and (I2.37|) are satisfied. 
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The paper is organized as follows. The next section deals with approximating problems. 
Theorem 12.41 is proved in Section H] and our argument relies on some a priori estimates on the 
approximate solutions and on monotonicity and compactness methods. Section [5] is devoted to 
prove Theorem 12.51 and uses a Moser type procedure. Finally, the proof of Theorem 12.71 is given 
in the last section. 

In our proofs, we use the well-known inequalities we are going to recall. As C M 3 is 
bounded and smooth, the Poincare inequality 

\\v\\v < -Mfj||Vu||# for every v G Vq (2.39) 

holds, and the space V is continuously embedded in L 6 (J7), i.e., 

IMIl 6 (Q) < -^olMlv for every v G V. (2.40) 

For the sake of completeness, we recall a related embedding result for parabolic spaces (see, e.g., 
|18l formula (3.2), p. 8]). For m > 1, we have 

L°°(0,T;L m (n))nL 2 (0,T;V ) Cl ?(ro) (Q) where g(m) := \ (m + 3) 
the embedding being continuous, i.e., 

IMlL<?(m)(Q) < ^,T,m(jMlL°°(0,T;L m (n)) + \\^ v \\l 2 {Q)^ 

for every v G L°°(0, T; L m {£l)) n L 2 (0, T; V ). (2.41) 
In particular, we observe that 

L°°(0, T; H) n L 2 (0, T; Jfc) C L 10 / 3 (Q) (2.42) 

L°°(0, T; L 3 (n)) n L 2 (0, T; Vo) C L A {Q) (2.43) 

and that the corresponding estimates (|2.4ip hold. 

The above inequalities are widely used in the sequel, as well as the elementary Young in- 
equality 

ab < 5a p + c 5jP V' Va,6>0 V<5>0 (2.44) 

where p,p' > 1 satisfy (1/p) + (1/p') = 1 and c<5 iP := {p')~ l {5p)~ p ' l p . 

We conclude this section by stating a general rule we use for constants, in order to avoid 
a boring notation. Throughout the paper, the symbol c stands for different constants which 
depend only on O, on the final time T, and on the constants and the norms of the functions 
involved in the assumptions of either our statements or our approximation. In particular, c is 
independent of the approximation parameter e we introduce in the next section. A notation like 
c<5 allows the constant to depend on the positive parameter 6, in addition. Hence, the meaning 
of c and c$ might change from line to line and even in the same chain of inequalities. On the 
contrary, we use different symbols (see, e.g., (|2.39[) ) to denote precise constants which we could 
refer to. By the way, several of such constants could be the same (like in (12.390 and in (I2.40|) ). 
since sharpness is not needed. 

3 Approximating problems 

This section contains a preliminar work in the direction of proving Theorem 12.41 and deals with 
a suitable appoximation of problem (|2.28|) - (|2.3ip . Namely, we replace the strong nonlinearities 
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that appear in equation (|2.28|) with smooth functions depending on the parameter e G (0,1). 
First of all, we see the logarithm In as a maximal monotone operator in R with domain (0, +oo). 
Precisely, for r£l, the set lnr is the singleton {lnr} if r > 0, while it is empty if r < 0. Then, 
we can consider the function Ln £ : R — > R defined as follows 

Ln £ r := er + ln £ r where ln e is the Yosida regularization of In. (3.1) 

We note that ln e is monotone and Lipschitz continuous with constant 1/e (see, e.g., |11| Prop. 2.6, 
p. 28]). Thus, e < Ln £ (r) < e + (1/e) for every r£l. The functions ln £ and Ln £ act on L 2 - 
spaces as well and Notation 12.31 is extended to them. Moreover, we replace F and G by new 
functions, still termed F and G, satisfying 

F' and G' are bounded (3.2) 

in addition to (|2.3p - (|2.4l) . Indeed, we can arbitrarily modify F and G outside [X*,X*] due 
to the last part of Theorem 12.41 Finally, we replace the possibly singular function (3 by a 
C°° function /3 e , in order to justify the chain rules we have to use. We proceed by extension, 
truncation, and regularization. For e > 0, we set 

n e :={xefl: dist(x,r) < e} and 1^. := {x e M 3 \ H : dist(x, V) < s}. 

As Q is of class C 2 , there exists £o £ (0, 1) such that, for every x G f2 £Q , there exists a unique 
point x G f2 £o satisfying 

x -\~ x 

x := — - — G r and x — x is orthogonal to T at x' (3-3) 

the correspondence x i— >■ x being a bi-Lipschitz diffeomorphism of class C 1 from Q' eo onto £l £Q . 
Then, we define £1 and the extension- by-reflection operator ~ : L°°(£l) — > L°°(ri) as follows 

Q := U £Q and for v G L°°(f2) and a. a. x G f2 we set 

v(x) = v(x) if x G and -u(x) = v(x) if x G r2 £Q . (3-4) 

Next, we extend further by reflection as well. We define Q and £ : L°°(Q) — > L°°(Q) as follows 

Q := O x (-T, 2T) and for v G L°°(Q) and a.a. (x, t) G Q we set 
(6t;)(aj ) t) = u(x,t) ift€(0,T), (£«)(x,t) = u(s,-t) ifiG(-T,0) 

and (£«) (x, t) = «(x, 2T - t) if t G (T, 2T) . (3.5) 

It is clear that the extension operator £ is linear and continuous. More precisely, we have 

sup ess Ev = sup ess v and infess £u = inf ess v for every v G L°°{Q). (3.6) 
q Q q Q 

Moreover, one can check that for every v G L°°(Q) we have 

Ev > a.e. in Q whenever v > a.e. in Q (3-7) 

HV^IUoo^) < M||Vt;|| L oo W) ifV^GL°°(Q) (3.8) 

WdtEvU^Q) < M\\dtv\\ L oo {Q) if 5^ G L°°(Q) (3.9) 

lip(£v) < Mlip-u if v is Lipschitz continuous (3.10) 

for some constant M, where lipt> is the Lipschitz constant of v. At this point, for e G (0,1) 
we define the operators $, (3 £ : Q x R — > R by using the extension operator £ and a truncation 
procedure as follows 

P(x,t,r) := (E/3(-,-,r))(x,t) (3.11) 
(3 E (x,t,r) := (3(x,t,r s ) where r £ := max{e, min{r, 1/e}} (3-12) 
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and observe that /3 £ is globally Lipschitz continuous. Indeed, by recalling (|2.5|) and setting for 
convenience 

L s :=lipP\ Qx[sm for 5 G (0,1) (3.13) 

we clearly see that ML £ is a Lipschitz constant for (3 £ . Moreover, as £ is linear and (|3.7|) 
holds, we infer that both P(x,t, ■ ) and P £ (x,t, ■ ) are nondecreasing on R for every (x,t) G Q. 
Furthermore, both (3 and /3 e vanish at r = 1. In particular, their values at every r € R have the 
sign of r — 1. 

Finally, we are ready to define a C°° approximation /J e : Q x R — > M of /3. We regularize /3 e 
by convolution and restrict the regularization we obtain to Q x R. Namely, we fix a nonnegative 
C G C°°(R 5 ) supported in the unit ball B of R 5 and normalized in L 1 (R 5 ). Then, by assuming 
£q < T and e G (0, eo) (such restrictions are not stressed in the sequel, but it is understood that 
they are satisfied), we recall A3. 13f) and set 

5 £ :=-^— and ( £ (x,t,r) := 5; 5 (((x,t,r)/5 £ ) for (x, t, r) G R 5 (3.14) 



P £ {x,t,r) := (0 £ * Ce)(x,t,r) = / £ (x - y,t - r,r - s) ( £ (y,T, s) dy dr ds 

JS E B 

= / (3 £ (x — 5 £ y,t — 5 £ T,r — S £ s) C(y,r, s) dy dr ds for (x, t, r) G Q x R. (3.15) 
Jb 

The reason of the above choice of 5 £ is that we would like to have 

\P e (x,t,r)-p e (x,t,r)\<Me for every (x, t, r) G Q x R (3.16) 

and some constant M. Actually, (|3.16p holds with the constant M that makes (|3.10p true, as 
we show at once. We have indeed 

\Pe(x,t,r) - P £ (x,t,r)\ = / (p £ (x - 5 £ y,t- 5 £ T,r - S £ s) - P £ (x,t,r)) ((y,T, s) dy dr ds 

Jb 

< I ML £ \{5 £ y, 8 £ t, 5 e s)\((y, r, s) dy dr ds < ML £ 5 £ < Me 
Jb 

since ML £ is a Lipschitz constant for (3 £ , as just observed. With a similar argument, we see that 

P e is Lipschitz continuous with lip/3 e < ML £ (3-17) 

since such a property holds for (3 e . In the sequel we use the following more precise facts 

sup \P E \ < _ sup \P\ and Hp P £ \ Qx [s,i/s] < ^pP\qx[8/2,i/8+8/2] 
Qx [5,1/5] Qx [5/2,1/5+5/2] 

for 5 6(0,1) and e<5/2. (3.18) 

Indeed, we have 5 £ < e < 5/2. Hence, if (x,t) G Q and 5 < r < 1/5, the values of p e in 
(|3.15p actually are values of (3 at points of the set Q x [5/2, 1/5 + 5/2], where P is bounded and 
Lipschitz continuous. Therefore, both the supremum and the Lipschitz constant are preserved 
by the convolution, since £ is normalized, and (|3.18p follow. Finally, we point out that 

P £ (x,t, • ) is nondecreasing on R for every (x,t) G Q (3.19) 

since such a property holds for (3 £ and £ is nonnegative. Moreover, we set for convenience 

P e (x,t, s) ds for a.a. (x,t) G Q and every r G R. (3.20) 
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Clearly, f3 E is convex with respect to the third variable. Furthermore, as P £ (x,t,r) and r — 1 
have the same sign as just observed, we see that (|3.16p implies 

Pe(x,t,r)>J^ ((3 £ (x,t,s) -Pe(x,t,s)) ds > -Me\r - 1\ for every (x,t,r) EQxR. (3.21) 
Finally, it is clear that the first of (|3.18j) implies the analogue for (3 e , namely 



sup \(3 £ \ < c s for every S G (0, 1) (3.22) 
)x [6,1/6] 



for some constant cs- 



At this point, we are ready to state the approximating problem, which consists in finding a 
triplet (t? £ ,X e ,£ £ ) having the proper regularity and satisfying 

$(0 e -G(X e ))-AI^0 e + £ e = 7r(0 e ) inL 2 (0,T;y *) and £ E = &(4) (3.23) 

d t X e - AX e + F'(X £ ) + G'{X £ ) 4 = a.e. in Q (3.24) 

4 = i? r a.e. on V x (0, T) (3.25) 

4(0) = 4 and X £ (0)=X . (3.26) 

The following result holds. 

Theorem 3.1. Let the assumptions of Theorem \2.J\ be fulfilled. Moreover, assume (|3.ip . (|3.2D . 
and ()3.15p . Then, problem (|3.23|) ( T3.26I) has a solution satisfying 

4 6 L 2 (0, T; V) n -ff 1 (0, T; Vq*) (3.27) 
X e € L 2 (0,T;H^) n IT^O, T; fl") (3.28) 

We avoid proving Theorem 13. II in order not to make the paper too long. Indeed, the a priori 
estimates in Section H] suggest how to proceed. Anyway, a rigorous proof could be done by 
regularizing ir and using, e.g., a Galerkin procedure. 

On the contrary, we prove some auxiliary results regarding the approximating nonlinear ities. 
The first formula we state is showed in [6l Lemma 6.1]. We repeat the short proof here, for 
convenience, and include it in the following proposition. 

Proposition 3.2. We have 

lnj 1 (s) = e s + es for every seR (3.29) 
In r 

- — < ln e r < lnr for every r > 1 (3.30) 

hn' £ (r) > 1 for every r < 1 and Ln^(r) > — for every r > 1 (3.31) 
4 < ln e r < t for every r G (3.32) 

where we have set I* := min{0, In 4} and i* := max{0, hit?*}. Moreover, we have 

2^ < K(0 < f- ewe n/ r e #1 ( 3 - 33 ) 

/or e small enough. 
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Proof. For r£l, let p £ (r) > be defined by the equation 

p e (r) +e\np £ (r) = r. (3.34) 
Take now any s € R. Then, r := e s + es satisfies p £ {r) = e s . On the other hand, we have 

ln e r= r -^l (3.35) 

by definition of Yosida regularization. We deduce that ln e r = s and (13.290 follows. To 
prove (|3,30p . we observe that for s > we have e £S > 1 + es and e s > 1. We infer that 

e (l+s)s > e s +£s > ^ 

If r > 1, applying this to s := ln e r (which is nonnegative since ln e 1 = 0), we obtain 

e (l+e) In. r > gln £ r + ^ r > gln £ r 

and (|3.30p follows from (|3.29)) by applying the logarithm. Next, we prove (I3.3ip . To this aim, 
we compute \n' e (r) from (|3.35|) and (|3.34p . We have 

;/ x l-p'Jr) 1 / Pe(r) \ 1 , 
ln' £ {r) = = - 1 - 7 ; = — . 3.36 

e £\ Pe{r)+e) p £ {r) + £ 

On the other hand, we observe that (|3.34j) and p £ {r) > 1 imply r > p £ (r) > 1. Hence, p £ (r) < 1 
for r < 1. We conclude that 

L4(r) =e + K(r) > e + = l \^f > 1 

for every r < 1. Assume now r > 1. Then, r + elnr > r, whence p e (r) < r. Accounting 
for (I3.36|) . we infer that 

i / \ 1 1 
In' (r) > > — 

eV ' ~ r + £ ~ 2r 

and the second of (|3.31|) follows. To prove (|3.32p . we observe that exp^* + £i* < and 
exp£* + £t ><&*. We deduce that exp4 + £t* < r < expt + £t for every r G [#*,#*], and 
(|3.32|) follows from (|3.29|) . Finally, we prove (|3.33|) . Owing to (|3.36|) and to the monotonicity 
of p £ , we see that r G implies that 

1 , #, s 1 
< In r) < 

On the other hand, it is clear that p £ {r') tends to r' as e — > for every r' > (see also, e.g., pTT| 
Prop. 2.6, p. 28]). Then, (I3.33|) immediately follows if e is small enough. □ 



Next, we deal with the analogue of (|2.9p for (5 £ . At the same time, we prove an inequality 
involving (5 £ . A notation like (|2.7p is extended to such functions. 

Proposition 3.3. We have 

\/3 £ , x (x,t,r)\ + \p E ,t(x,t,r)\ <c(l + \/3 e (x,t,r)\) and \%, t (x,t,s)\ < c{%(r) + |r| + l) (3.37) 
for every (x, i,r)eQxR, some constant c, and £ small enough. 
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Proof. As far as the first inequality is concerned, we deal, e.g., with /3 E t, since the argument for 
the space derivatives is quite similar. We first prove that 



\d t %{x,t,r)\<MM p (l + \%{x,t,r)\) 



(3.38) 



for every (x, t) 6 Q and r E E, where M and Mp are the constants satisfying (|3.8|) - (|3.10p and 
(|2.9|) . respectively. Assume first r G [e, 1/e] . Then, (|3.38p coincides with the analogue for (3 due 
to (|3.12p . and this easily follows from ()3.9p . (|3.1ip . and ()2.9p . Assume now r < e. Then, we have 

\dtfi e (x,t,r)\ = \d${x,t,e)\<MM l3 (l + 0(x,t,e)\) = MM p (l + \%(x, t, r)|) 

just using the above with r = e. As the argument is similar if r > 1 /e, (|3.38p is established. In 
order to prove the first of (13.37p . we notice that 

± P E ,t = ±d t {% * Ce) = ±($&) * Ce < c(l + \0 e \) * Q = C+ C%\ * Q 

with c := MM p, since (I3.38H holds and the convolution with the nonnegative normalized kernel 
Q preserves order and constants. Therefore, we have to bound the last convolution with the 
right-hand side of the inequality we want to prove. Assume first (x, t) € Q and r > 1 + e. Then, 
(3 £ (y, r, s) > for (y, r) € Q and \s — r\ < 5 £ (since S £ < e), and we have 

($ e \*Ce)(x,t,r) = (P £ *( £ ){x,t,r)=(3 £ {x,t,r) = \(3 £ {x,t,r)\. 

If r < 1 — e the argument is similar. Finally, if \r — 1| < e, by assuming e < 1/4, we have 
5 £ < 1/4, r — 5 £ > max{e, 1/2}, and r + 5 £ < min{l/e, 3/2}, whence 



\(\P £ \*( £ )(x,t,r)\ 



* Ce)(x,t,r)\ 



— _ SU P 

Qx[r-S e ,r+S e ] 



< 



sup 

)x [1/2,3/2] 



since C £ is normalized in L 1 (M 5 ), and the first of (]3.37p clearly follows. Then, we easily derive 
the second inequality. As (3 £ is monotone with respect to the third variable, we have 



\%,t{x,t,r)\ 



/3 £ j(x,t,s)ds 



ds 


< c 


I 









'l + \(3 £ (x,t,s)\)ds 



<c\r-l\ + c J (P e (x,t,s) - P e (x,t,l)) ds + c 
<c(%(x,t,r) + \r\ + l) 
since (3 £ (x,t, 1) is bounded due to (13. 18H with 5 = 1/2. 



\Pe(x,t,l)\ds 



□ 



The last preliminary remarks we make regard relationships between the approximating non- 
linearities and the boundary datum $r- In the next section, we need to consider a known smooth 
function that coincides with $r on the boundary. Thus, a natural choice is the harmonic exten- 
sion of i?r- Precisely, we define fi-K : Q — > M by the conditions 

ifj{(t)e7, tfw(i)|r = #r(*), and M K (t) = in O, for a. a. t G (0, T). (3.39) 

Then, assumption (I2.20p and the general theory of harmonic functions (in particular, the maxi- 
mum principle) ensure that £ L 2 (0, T; y)nif 1 (0, T; H) and that the following estimates hold 

\\$x\\L 2 (0,T;V)nm(0,T;H) < C \\^T \\L 2 (Q,T;H 1 / 2 (r))nH 1 (Q,T;H- 1 / 2 (r)) 

and ^ < tin < $* a.e. in Q. (3.40) 
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Proposition 3.4. We have 

l|Lne#rK||L»(Q)n£ a (0,T;V) < c (3.41) 
H/5e(^)|lL°°(Q)nL 2 (0,T ; y)n/f 1 (0,T;/f) < c ( 3 -42) 

/or some constant c and e small enough. 

Proof. Estimate (|3.40p and (13.321) imply that i* < ln e t?^ < I* a.e. in Q, where the notation of 
Proposition 13.21 has been used. Owing to (|3.33|) as well, (|3.41|) immediately follows. To prove 
the L°°-bound of 4332]), it suffices to recall p72|) . (I2TT3D . and the first of (ETTKl) . Finally, we 
prove the estimate regarding the time derivative (the argument for space derivatives is similar). 
We have 

11^(^)11^(0) < \\Pe,t&x)Wwm + \\P'e(#x)dtM\LHQ) 

<c\\P e ,t(fix)\\L°°(Q)+ SUP \P' e \ Wrf L? (Q) < c 
Qx[i,i)«] 

due to (I3.37|) . the L°°-estimate just proved, the second of A3. 18H . and ()3.40p . provided that e is 
small enough. □ 



4 Existence 

In this section, we prove Theorem 12.41 More precisely, we first check the last part of the 
statement, i.e., we show that every solution to problem (|2.28p - (l2.3ip satisfies the bounds p. 321) . 
Then, we consider the approximating problem (|3.23p - (|3.26p taking the further assumption (|3.2p 
into account, and perform a number of a priori estimates on its solution (i9 e , X e , £ e ). Finally, we 
let e tend to zero by using monotonicity and compactness methods. This leads to the proof of 
Theorem 12.41 under the additional condition (|3.2|) . However, it is clear that such a procedure 
proves Theorem 12.41 in the general case. Indeed, the triplet X, £) we find solves anyone of 
the problems obtained by replacing the original functions F and G by other ones that coincides 
with the given F and G on [X*,X*], due to 07321) . 

Hence, we assume (13. 2p throghout the present section. As such a condition is added to our 
assumptions, we allow the values of the different constants c to depend even on the Lipschitz 
constants of F and G, following the general rule explained at the end of Section [2j Moreover, 
it is understood that e belongs to (0, 1) and is even smaller as in the previous section (see, e.g., 
Propositions 13.31 and 13. 4|) . Finally, 5 is a positive parameter, say 6 £ (0, 1), whose value is chosen 
according to our convenience. 

The maximum principle argument. We prove that every solution (-#, X,£) to problem 
(pT28]) - (p73T]) satisfies (pQ2|) . Fix a Lipschitz continuous function H : R -> E of class C 1 
such that 

H(r)=0 ifre[X*,X*] and H'(r)>0 if r £ [X*, X*]. (4.1) 

As X G L 2 (0, T; V), it turns out that v := H(X) is an admissible test function for (I2.29p . Then, 
we multiply (|2.29|) by H (X) and integrate over Qt, where t € (0, T] is arbitray. After integrating 
by parts, we obtain 

f d t XH(X)+ [ VX-VH(X)+ / (-F'(X) + G'(X)tf) ff(X) = 0. (4.2) 
JQt JQt JQt 
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If H denotes the primitive of H vanishing at X*, we have 

[ d t XH(X)= [ H(X(t))- [ H(Xq)= ( H(X(t))>0 
JQt Jn Jn Jn 

due to (|2.22|) and (|4.1|) , The next integral of (|4.2|) is nonnegative, obviously, and the last one 
has the same property. Indeed, its integrand is nonnegative since i? > and the three functions 
F', G', and H have the same sign, due to (|2.4|) and (|4.1|) . Therefore, all the integrals of (|4.2|) 
vanish identically. In particular, we deduce that H(X(t)) = a.e. in Q, for every t 6 [0, T], and 
this clearly implies (|2.32D . 

Now, we start proving estimates on the solution ($ £ ,X £ ,£ £ ) to problem (|3.23|) - (|3.26|) . as 
mentioned at the beginning of the present section. When we shortly say that we test an equation 
by some function v, we mean that we test such an equation at time s by v(s), we integrate first 
over O and then over (0, t) with respect to s, where t € (0, T) is arbitrary, and we integrate by 
parts, if necessary. 

First a priori estimate. We test (|3.23p by v := $ e — dji + <5(Ln £ -d £ — Ln £ 7%c). At the same 
time, we test (|3,24p by dtX £ . Then, we sum the equalities we get to each other and note that 
two terms involving G cancel. After rearranging a little and adding the same integral to both 
sides for convenience, we obtain 



\ [ \Mt)-#x(t)\ 2 + I dtG(X e )#M+ [ VLn £ tf £ -W, 
z Jn JQt JQt 

JQt 

+ 5 ( dt# e Ln £ $ £ + 5 [ | V(Ln £ e - Ln £ ij^)\ 2 
JQt JQt 

+ 5 [ - /3 £ 0%t))(Ln e #e ~ Ln £ #x) 

JQt 

+ [ \dtX £ \ 2 + l [ |VX £ (t)| 2 + i / \X £ (t)\ 2 
JQt 1 Jn * Jo. 

|^o-^(0)| 



/ dtfnWe - + ! V Ln £ $ £ ■ 
JQt JQt 

JQt 

+ 5 [ d t $ E Ln £ 0m + S [ G'(X £ ) d t X £ (Ln £ e - Ln £ 
JQt JQt 

-5 [ V Ln £ tin ■ V (Ln £ 4 - Ln £ K ) - 5 / & (^) (Ln E e - Ln £ 
•/<2t ■/ Qt 

+ J / |VX | 2 + 5 / vr(4)(Ln £ 4-Ln £ ^)- / F(X £ )d t X £ + i / |X £ (t)| 2 . (4.3) 
1 Jn JQ t JQt 1 Jn 

Now, we observe with the help of monotonicity that all terms on the left-hand side are nonneg- 
ative but two of them. We deal with the first one that needs some treatment. We have 

d t G(X £ )$x= [ G(X £ (t))^(t) - [ G(Xo)4k(0)- f G(X £ )d t tfx 
Jn Jn JQ t 

> tf* / G(X £ (t)) -c-\l |G(X £ )| 2 -\! \d t $x\ 2 f G(X £ (t)) -c [ |X £ | 2 - c. 
Jn 1 JQ t 1 JQ t Jn JQt 



The second one is the following 



5 I d t # E Ln £ tf £ = 5 / H e (Mt)) - 5 / L £ ($ ) > 5 / £ e (0 e (t)) " c 
Jn Jn Jn 



where we have set 



/r g /•)■ 

Ln e s ds = - (r 2 - 1) + y ln £ s ds for r € R. 

Note that £ £ is convex and bounded from below uniformly with respect to e since ln E 1 
Now, we consider the right-hand side and deal with the non-trivial terms of it. We have 

- / dfdu@ s - M ) < \ f \d t ^\ 2 + [ |4 " $K\ 2 < f % ~ $K\ 2 + C 
JQt 4 JQt JQt JQt 

Next, we consider 

( VLn £ 4-V^= ( V(Ln e i? e -Ln e ^)-V^+ / VLn^M-Ww 
JQt JQt JQt 

< 5 - [ |V(Ln e 4-Ln £ ^)| 2 + ^ / |V^| 2 + c < 6 - [ |V(Ln £ e - Ln £ 4k)| 2 + c s 
8 JQt d JQt 8 JQt 

Moreover, by (|3.42p and (|2.1ip it is easily seen that 

WnW* - 4k) + / vr(^)(4 - K ) 
<c/ (l + |4| + ko|) |4-4k| <c / |4-4x| 2 + c 

JQt J Qt 

We deal with the next integral of (|4.3|) as follows (cf. (|3,41|) ) 

5 [ d t ti £ Ln £ #x = S [ dt(0 e - K ) Ln £ ^ + 5 / Ln £ 

JQt JQt JQt 

= 5 [ (A(t)-#x(t))Ln £ #x(t)-6 [ (^ -^(0))Ln £ ^(0) 
Jn Jn 

-6 f (4 - #H)dt Ln £ 4k + <5 / Ln £ TC 

JQt JQt 

/ [0 6 (t)-0 M (*)| 2 + / |4-4kI 2 + c. 

2 J^ JQt 
Subsequently, by the Poincare inequality (|2.39p we have 

5 [ G'(X £ )d t X £ (Ln £ 4 - Ln £ fa) 
JQt 

~8lP J 1 Lll£ 4 " Lll£ ^ + 2 * M " / l 6 ^ 9 ** 

< | / |V(Ln £ 4 - Ln £ 4k) | 2 + 25M^ sup |G"| 2 / |$X, 
8 JQt JQt 

The next integral is easy. We have indeed 

-S f VLn £ ^. V(Ln £ 4 - Ln £ 4k) < | / | V(Ln £ 4 - Ln £ 4k)| 2 + c 
JQt 8 JQt 



2 



2 

£ | • 
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while the other two terms involving Ln £ are estimated owing to the Poincare inequality once 
more this way 

S [ /? £ (4k)(Lu £ 4 ~ L% 4k) + 6 [ 7r(4)(Ln £ 4 - Ln £ 4k) 
<-zLs I \Im e # e --Ln e # M \ 2 + 26M% f |tt(4) - &(4k)| 2 

|V(Ln £ 4-Ln e ^)| 2 + c / |tt(4)| 2 + c/ |&(4k)| 2 
8 JQt JQ* JQt 

<|/ |V(Ln £ 4-Ln e ^)| 2 + c / |4 - 4k| 2 + c. 
8 JQt JQt 

Next, we treat the second last integral of (|4.3p . We have 

-/ F{x £ )d t x £ <- f \d t x £ \ 2 + [ \F{X £ )\ 2 <-f \d t X £ \ 2 + c[ \X £ \ 2 + c. 
JQt 4 JQt JQt 4 JQt JQt 

Finally, we deal with the last term as follows 

\ [ \Mt)\ 2 = \ [ \M 2 + [ *ed t X £ <] [ \d t X £ \ 2 + [ \X £ \ 2 + c. 
z Jn z Jn JQt 4 JQt JQt 

Hence, we can recall (]4.3f) and all the estimates we have derived. Forgetting some nonnegative 
terms on the left-hand side we have 

' [ |4(*)-4k(*)| 2 + [ VLn £ 4"V4 + ^ / |V(Ln £ 4-Ln £ ^)| 2 
Jn JQt 1 JQt 



2 

+ Q-25M 2 snp|G'| 2 ) I \d t X £ \ 2 + \ f \VX £ (t)\ 2 + \ f |X £ (i)f 



JQt JQt 



Therefore, we can choose 5 small enough and apply the Gronwall lemma. By the Poincare 
inequality, we conclude that 

114 - &m\\l°°(o,T;H) + ll( Lr 4(^)) 1/2 We\\l?(0,T;H) + ||Ln £ 4 -Ln £ 4k||i,2(0,T;V) 
+ ll^e||L 2 (0,T;fZ") + ll X e||L°°(0,T;V) < c - 

Then, thanks to (|3.40|) and (|3,41|) we obtain the basic a priori estimate 

H4||L~(o,T;tf) + ll(Ln £ (4)) 1/2 V4II 

+ l|Ln £ 4||L 2 (0,T ; y) + \\ X e\\L^(0,T;V)nm(0,T;H) < c - ( 4 -5) 

First consequence. As all the terms of (|3.24p but the Laplacian are bounded in L 2 (0,T; H) 
due to (14. 5p and (13.21) . we immediately derive that the Laplacian is bounded as well, whence 



\\ X s\\l 2 (0,T;W) < c ( 4 - 6 ) 

by the homogeneous Neumann boundary condition satisfied by X £ and the general theory of 
elliptic equations. 
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Second consequence. We set Qf := {(x,t) G Q : $ £ (x,t) > 1} and Q~ := {(x,t) G Q : 
$£{ x ,t) < 1}. Moreover, if t> is a real function on Q, v 1 * 1 := max{±?;,0} are its positive and 
negative parts. Then, inequalities (|3.3ip yield 



|2 



/ L<(4)|V4| 2 > / L<(4)|W,| 2 > / |W £ 
JQ7 JQ7 

whence (|4.5|) implies that 

|Vtf £ | 2 <c and / |W £ 1/2 | 2 <c. (4.7) 



Now, we set n~(t) := {x G Q : e (t) < 1} and Sl+(t) :={x£0: e (i) > 1} for a.a. t G (0,T) 
and use (|4.7f) in order to estimate Vt? e in a suitable norm. Accounting for the first (|4.7f) . we 
see that 

[ T \m(t)\\l* /S(Q7{t)) dt < c [ T \m(t)\\ 2 LHQ7m dt = c[ < c . (4.8) 

On the other hand, we have W e = 2$l /2 \7$l /2 a.e. in Q+. Therefore, using the Holder inequality, 
we obtain for a.a. t G (0, T) 

l|W«(t)|| £ 4,.(n+(t,) ^ 2 ll^ 1/2 (*)llL*(n+(t))ll^ e 1/2 («)llL a(n + W) 

<2||4(t)|li / 2 2 (f , ) ||v^/ 2( i)|| i2( ^ {t)) 

and squaring, integrating over (0, T), and owing to (14. 5p and to the second (|4.7|) . we derive that 

[ T \\VMt)\\% /Hnt{t)) dt < 4||4||loo(o,t ; //) f |V^/ 2 | 2 < c. (4.9) 
Finally, we observe that for a.a. t G (0, T) we have 

HV4(t)|| 4 / 4 3 /3(0) = \m(t)C HQm) + HV4(t)|| 4 / 4 3 /3(n , W) whence 

llv^WII 2 ^^ < c(||V4(i)|| 2 4/3 (nrW) + IIV4(t)lli 4/3(Q + (i)) ) 

so that, using (14. 8p and (14. 9p . we conclude that 

II V ^IIl2(0,T;L4/3 (C)) < C. (4.10) 

Second a priori estimate. We observe that 

(dfd e ) Pe(A) = dSei'&e) ~ %, t ($e) ^ V&(4) = /M^) + W e • 

Therefore, if we test (13.231) by f3 £ {$e) — f3e($Ji), rearrange, and add the same quantity to both 
sides in order to take advantage of (|3.2ip . we obtain 

(&(i,4(i))+M £ |4(i)-l|)+ / VLn £ ^-^(4)V4+ / - fe^x)) 2 

JQt JQt 

%(0, tfo) + / %, t (#e) + [ d t V E /3 e (0M) + Me [ \0 e (t) - 1| 
JQt JQt Jtt 

- [ VLn £ tf £ . (3 £ , x (4) + / V Ln e 4 • V/3 £ (tf M ) 

+ / (d t G(X £ ) - p e (&n) + *{#*)) (&(0 e ) - A^))- (4.11) 
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All the terms on the left-hand side are nonnegative, and we estimate each integral on the right- 
hand side, separately. The first term is bounded due to (|2.2ip and (|3.22p . The next one is 
treated by accounting for the second of (|3.37p this way 

/ %M<c( + |4l + i) 

J Qt JQt 

< C f %(#e)+C<c[ (%{$ £ )+Me%-l\)+C 

JQt JQt 

since i? e has already been estimated in L°°(0, T; H) by (|4.5p . For that reason, the fourth integral 
of (14.111) is bounded as well and we deal with the third one. We integrate it by parts and have 

/ ^ eJ 8 e (lJ M )= / «? e (t)ft(i,rfjt(i))- / ^(0,tf H (0))- / tiedtPetfx) 

JQt Jn Jn JQ t 

and one immediately sees that the last right-hand side is bounded, due to (|4.5|) and (|3.42|) . For 
the same reason, the next integral of (|4.1ip that involves is bounded as well. Now, we treat 
the term continainig f3 s>x - By using the first of (|3.37p and accounting for (|4.5p and (|3.42p once 
more, we easily have 

-/ VLn e tf £ -/Mtf e )<5 / {\Pe(A)\ + if + c s < 5 [ (&(0 e ) ~ &(*k)) 2 + c 5 . 
JQt JQt JQt 

Finally, owing to (pHj) . ((332]), and (i2~TT|) . we obtain 

/ (d t G(X E ) - f3 e (# x ) + 7T(4)) We{A) ~ Pe^Dl)) 
JQt 

<S [ \f3M) ~ Pe($x)\ 2 + C8 [ (\d t X e \ 2 + \$ e \ 2 + l) 
JQt JQt 

and the last integral is bounded due to (14. 5p . Therefore, if we collect (I4.1ip and the inequalities 
we have proved, choose 5 small enough, and apply the Gronwall lemma, we derive that 

\\Pe(# e )\\v(p,T;H)<C (4.12) 



Consequence. Estimates (|4.5p and (|4.12|) and our assumptions (|3.2|) and (|2.11|) on G and ir 
ensure that 

II<W|l2(0,T;V *) < c ( 4 - 13 ) 

just by comparison in (13.231) . 

Lemma 4.1. Assume z, z n € L 2 (Q), z > a.e. in Q, and z n — > z a.e. in Q. Moreover, let {e n } 
be a positive real sequence converging to 0. Then, {Pe„( z n)} converges to /3(z) a.e. in Q. 

Proof. It suffices to show that, for every 5 £ (0, 1), we have 

0e„{z n ) -> P(z) almost uniformly in Q s := {(x,t) 6 Q : S < z(x,t) < 1/6}. 

Thus, fix 5 € (0, 1) and rj > 0. We have to show that a subset C Q s exists such that 
\Q\Qt\ < V and (3 £n (z n ) — > j3{z) uniformly in Q 5 , where | • | stands for the Lebesgue measure 
in M 4 . By the Severini-Egorov theorem, we find C Q s such that \Q \ Q*L\ < rj and z n — > z 
uniformly in Q^, and we can prove that (3 £n (z n ) — > (3{z) uniformly in Qf r Fix n such that 

S S 2 

£n < 7: and - < z n < - in Q 5 for every n > n. 
v v A ' 
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On the other hand, we have 

\\[3e n {Zn)-P{z)\\ L oo {QSn) < \\PeAZn)-peM)\\L°°(QS) + \\Pe n (z)-P(z)\\L°°(QS)- (4-14) 

Assume now n > n. Then, e n < 5, whence e n < z < l/e n . Thus, j3(z) = $ £n (z). We infer that 
the last term of (|4.14p is < Me n by (|3.16p . On the other hand, as e n < 5/2, we can use the 
second of (|3.18p . Therefore, we conclude that 

\\Pe n {Zn) ~ M\\l°°{Q*,) ^ °s\\ z n ~ A\l^{Q^) + M ?n 

and deduce that (3 £n {z n ) converges to (3(z) uniformly in Q^. □ 

Conclusion of the proof. The estimates (f4~5|) . (jOj) . KlO\i . (fiT2|) . and P~T3|) proved in the 
previous steps and classical weak and weak star compacness results ensure that suitable limit 
functions exist in order that the following convergences hold (at least for a subsequence) 

weakly star in L°° (0, T; H) n il 1 (0, T; Vq) (4.15) 

weakly in L 2 (0,T;V) (4.16) 

weakly in L 2 (0, T; W) n H l (0, T; H) (4.17) 

weakly in L 2 (0, T; L 4 / 3 (Q) ) (4.18) 

weakly in L 2 (0,T; H). (4.19) 

Note that P~T5]) and P~Tgj) imply that 4 converges to •& weakly in L 2 (0, T; W 1 ' 4 / 3 ^)). Now, 
we observe that the Sobolev exponent (4/3)* of W 1A/3 (Q) is 12/5 > 2. Hence, H /1 ' 4 / 3 (0) is 
compactly embedded in H. On the other hand, even W is compactly embedded in H. Therefore, 
by applying [27, Thm. 5.1, p. 58] and possibly taking another subsequence, we derive that 

4 -► ■& and X £ -> X strongly in L 2 (0, T; F) and a.e. in Q. (4.20) 

This allows us to identify all the limits of the nonlinear terms. As far as the logarithm is 
concerned, we note that (|4.15|) and (|4.16|) imply that ln e $ converges to £ weakly in L 2 (0,T;H). 
Hence, we can conclude that $ > and I = hit? a.e. in Q (see, e.g., [11} Prop. 2.5, p. 27] for a 
similar result). From $ > a.e. in Q and (|4.2U|) for we see that we can apply Lemma [4. II and 
infer that ^ = /3(i?) a.e. in Q. Finally, the limits of the remaining nonlinear terms (i.e., those 
related to G, F' , G', and ir) can be identified by using the convergences a.e. given by (|4.20p and 
accounting for our assumptions (12. 3D and p. lip . This concludes the proof of Theorem 12.41 

5 Boundedness 

In this section, we prove Theorem 12.51 by estimating the L p -norm of t? (or of a suitable function 
of it) with a constant independent of p by using a Moser type technique. As usual, if z is either 
a function or a real number, the symbol z + denotes its positive part. Moreover, it is understood 
that n is a positive integer and 5 is a positive parameter, say 5 & (0, 1). Furthermore, we set for 
convenience 

u* --lnti* (5.1) 

Finally, we assume q < 4 (see (|2.33p ) without loss of generality. 

In the sequel, we perform two a priori estimates. In each of them, the use of the chain rule 
for time derivatives has to be justified, and the trouble is a lack of regularity for drf, which is 



4 

Ln e 4 i 

X £ ► X 
W e -» W 
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not known to belong to L 2 (0,T; H). The first lemma we prove overcomes such a difficulty and 
uses just the regularity we actually know for namely 

i? G L°°(0, T; H) fl -ff 1 (0, T; Vq), •& > a.e. in Q, and In 7? G L 2 (0, T; V) (5.2) 

the regularity for the time derivative being a consequence of Theorem 12.41 (see Remark 12. 6p . 
Actually, the last of (|5.2[) does not play a special role in the lemma, in which a general continuous 
increasing function </>' : (0, +oo) — > ]R is considered. 

Lemma 5.1. Assume i9 G L°°(0,T;H) n -H' 1 (0, T; Vq*) and 7? > a.e. m Q. Moreover, let 
4> ■ (0, +oo) — > R 6e a convex function of class C 1 and assume that 0'(i?) G L 2 (0, T; Vb). Then, 
z/<3? : R — ► (— oo,+oo] denotes the extension 

$(r) := 0(r) if r > 0, <J>(0) := lim 0(r), and <&(r) := +oo if r < (5.3) 

j — »o+ 

i/ien £/ie function t i— ► <£($(£)) is absolutely continuous on [0,T] and we have 

[ (d t tf{s),(j)'(${s)))ds = [ $(#(*))- / $(tf(0)) /or every i G [0, T]. (5.4) 
Jo Jn Jc 

Proof. We first observe that $ is convex, proper, and lower semicontinuous in R. In addition, 
we notice that 

4>{u) G <9<3?(u) a.e. in O if u £ H, u > a.e. in 0, and 0'(tt) G if (5.5) 

and the conjugate function of <1> satisfies 9<$* = (<9$) -1 . For a. a. £ G (0, T) we set for 
convenience v(t) := </>'(i?(i)) and observe that both and v(t) lie in if, and i?(t) > a.e. in Jl. 
Thus, it turns out that v(t) G <9<l?(t?(t)) a.e. in £2 by (|5.5p . and consequently 

tf(t) G d$*(v(f)) a.e. in (5.6) 

Moreover, defining the functionals J : H — > (— oo, +oo] and Jo : Vb — > (— oo, +oo] as follows 

J(v):=[&*(v) if **(«) G L 1 ^), J(w):=+c» otherwise, and J := J|y (5.7) 

we note that (cf., e.g., [U Prop. 2.8, p 71]) J is convex, proper, and l.s.c, and its subdifferential 
operator dJ : H — > 2 H is exactly induced by via the almost everywhere in Q inclusion. 
Then, for a.a. t G (0,T) I^M entails $(t) G dJ(v(t)), that is, 

(0(t),v(t) + J(v(t)) < J(w) Vweff. 

Therefore, as v(t) G Vb (and v(t) G D(dJ) C D(J)), the functional Jo is proper and 

(0(t), u(t) - + J «t)) < JoH Vw G V 

whence the inclusion $(i) G dJo(v(t)) holds for the subdifferential <9Jo : Vq — > 2^o* as well. Then, 
introducing the conjugate functionals and the subdifferentials 

J* :H -» (-oo, +oo], J * : V * -» (-oo, +oo] and dJ* : H ^ 2 H , dJ * : V * -► 2 y ° 

and recalling that <9J* = (5J)" 1 and <9J * = (SJo)" 1 , we observe that G 9J*(i9(i)) and 
i>(i) G <9J(*(i?(t)) for a.a. i G (0,T). Therefore, being understood that i? denotes the Vq*— valued 
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continuous representative, from the latter we conclude that (see, e.g., Lemma 3.3, p. 73] for 
a similar result) Jq(i?) is absolutely continuous in [0, T] and 

\d t #(s),ct)'(#(s)))ds= I (d^(s),v(s))ds = J^(t))-Jo(m) Vte[0,T]. (5.8) 
o Jo 

Moreover, the same representative t? is H— valued and weakly continuous from [0, T] to H. Now, 
as $** = <3? we remind that for u, w G H one has 

J*( u ) ■= ( $(u) if G J* (u) := +oo otherwise (5.9) 

Jn 

w G dJ*(u) if and only if w G <9<I>(u) a.e. in 0. (5.10) 
We also claim that 

J*(u) = J*(u) if ueff. (5.11) 

For n G -ff we have indeed 

Jq(u) = sup{(u,t(;) — Jo('" ; )} < sup {(ii,u;) — J(w)} = J*(u). 

On the other hand, in view of |17l Lemma 2.3] (or, also for related results, [2J Lemma 2.4 and 
Section 2]), it turns out that for all w G D(J) there exixts a sequence {w n } C Vo such that 
w n — > w in H and Jo(w n ) = J(w n ) — > J(w) as n -> oo, whence 

(u,w) - J(io) = lim {(«,«;„) - J(w n )} = lim {{u,w n ) - J (w n )} < J *(u) 

n— »oo n— »oo 

and consequently J*(u) < Jq(u) as the inequality (u,w) — J(w) < Jq(u) holds for all w G H. 
Then, (|5.1ip is proved. By combining (|5.9[) and (|5.1ip we conclude that 

j*mt)) = j*(m)= [ nm) vteM. 

Jn 

This yields the assertion of the lemma. □ 



Remark 5.2. We can replace <3> by (ft in the right-hand side of (15. 4ft in a number of cases. For 
instance, if we know that (the continuous representative of) i? satisfies t?(i) > a.e. in 17 for 
every t G [0, T], then <J>(i?(i)) = (f)(${t)) a.e. in Q for every t as well. Similarly, we have the same 
conclusion, independently of the strict positivity of whenever <£(0) = +oo. 

Lemma 5.3. Set 

<j) n {r) := ^ e 2min{„,(lns-«*)+} _ ^ ds ^ f £ ( 0)+oo ). (5 12 ) 

Then, 'positive constants a* and C* exist such that 

Mr) > e 3min{n ' (lnr - lt * )+} - C* (5.13) 
for every r G (0, +oo) and any positive integer n. 

Proof. Assume first $* < r < r d*e n . Then, we have 
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for some C* > 0, whence a* := i?*/6 works in (|5.13|) for this case. Moreover, we can assume 
C* > a*, so that (15.131) holds even for r £ (0,??*), since 4> n (r) = for such values of r. Finally, 
if r > #*e n , we have r > r' := $*e n and we already now that (|5.13p holds with r = r'. We 
deduce that 

Mr) > Mr') > a* e 3min K( lnr '-«*) + } - C* = a* e 3n - C* = a* e 3min{n,(inr-«*)+} _ c * 
This concludes the proof. □ 
Lemma 5.4. Assume p £ [1, +oo) and set 

i]) n (r):= i (min{n, (Ins - u*) + } 2p ~ l ) ds for r £ (0, +oo). (5.14) 



Then, we have 



^n(r) > — min{n, (lnr — u*) + } 2p for every r £ (0,+oo). (5.15) 
2p 



Mr)= I ey{y-u*fP- 1 dy> I (y - u*) 2p ~ l dy = ^(lnr - u*) 2p 



Proof. If -d* < r < $*e n , we have 

flnr /-lnr 

(„ _ „*^2P-1 di, = 

2p 

and ()5.15|) follow. If instead, r > i?*e n , we observe that ij)n( r ) > '4 ) n{ r ')i where r' := i?*e™ '. On 
the other hand, we already now that (|5.15p holds with r = r'. Hence, we easily conclude that 
the desired inequality is true even in this case. Finally, if r < i?*, we have ip n (r) = and (I5.15D 
trivially holds. □ 

Now, we start estimating. 
First a priori estimate. We set 

u:=ln$, w n := min{n, (u — u*) + }, and v n := e 2w " — 1. (5.16) 

We want to use v n as a test function in (|2.28[) and apply Lemma 15.11 with <f> = <j) n given 
by (|5.12p . To this aim, we note that u £ L 2 (0, T; V) and that v n = (j> n (u), where <f> n is a Lip- 
schitz continuous function. Hence, v n £ L 2 (0,T;V). Moreover, v n vanishes on the boundary 
since ln^r < u* by (|2.20p . Therefore, v n £ L 2 (0,T;V ). Furthermore, 4> 

n is a C convex 

function on (0, +oo) and v n = 4>' n {'&)- Hence, we are allowed both to test (I2.28P by v n and to 
apply Lemma [5.11 Note that </> n (i?o) = since #o < #* by (|2.2ip . Then, by extending 4> n (r) 
with value for r = (cf. (|5.3p ) and setting 

/ : = -f3{r ) + tt(#) + d t G(X) (5.17) 

from (15. 4p we obtain 

I M$(t))+ I Vu-Vv n + [ (P(#)-P(r))v n = I fv n . (5.18) 
Jn JQt JQt JQt 

We treat each integral, separately. For the first one, we have 

M^(t)) > "* / e 3M, " {t) -C* = a4e w "W\\ 3 L3m - c 



27 



thanks to Lemma |5 .31 The next term of (|5.18|) is treated this way 

( Vu-Vv n = [ Vw n -Vv n = 2 [ e 2Wn \Vw n \ 2 = 2 [ \Ve Wnl2 
JQt JQt JQt JQt 



and the last term on the left-hand side is nonnegative too. Indeed, v n is nonnegative and 
> where v n > since (5 is monotone. Thus, let us consider the right-hand side. We 

first notice that / G L 2 (Q). Therefore, using the Holder, Sobolev, and Poincare inequalities 
(see (jZHl) and (|ZMJ) ), we have 



fv n < 



J O.t 



<II/IUq)+ ril/(^)llL 2 (m||e^ (s) || L 3(n)l|e^ (s) 



Li(Q) -I- / \\J{S)\\L2(n)\\e \\L3{n)\\e Wi^n) 
Jo 

t pt 



<c + 5j o \\e^ s X Hn) ds + c 5 ||/(*)||2r||e«-W||i, (n) cfa 
< 28 f\e w ^ - 1||| 6(Q) ds + c 5 J*\\f(s)\\ 2 H \\e w ^\\ 2 L3{n) ds + c 8 



At this point, we collect all the above estimates, choose 5 small enough, and apply the Gronwall 
lemma (see, e.g., [TO Lemma A. 4, p. 156]) noting that & ^ 1 (0,T) since / G L 2 (Q). 

We obtain 

11^1^0.^3^+ / |Ve-"| 2 <c. (5.20) 



Q 



Consequence. From (I5.20p and ()2.43p . we deduce that exp(w; n ) is bounded in L 4 (Q). Hence, 
we can let n tend to infinity and infer that exp((u — u*) + ) G L 4 (Q). Now, we observe that 

e («-n*)+ = e -u* e u = fijft* where § > r _ 

As •& is positive, we conclude that 

G L 4 (Q). (5.21) 

Now, we rewrite equation (|2.29|) in the form 

d t X — AX = -F'{X) - G'(X) i? 

and observe that the right-hand side belongs to L 4 (Q) due to (|5.2ip and the Lipschitz conti- 
nuity of F and G. By the general theory for parabolic equations, we infer that dtX G L 4 (Q), 
whence also 

\\dtG(X)\\ Lq{Q) < c\\d t G(X)\\ LHQ) < +oo (5.22) 

since q < 4, as we have assumed at the beginning of this section. On the other hand, we can 
estimate the right-hand side of (|2.28p in a better way owing to (|4.12|) . Indeed, recalling (|2.33l) , 
we conclude that 

lkWlli9(Q) < c II^IIl4(q) + IKo|k«(Q) < +oo. (5.23) 
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The Moser type procedure. Our aim is to prove an iterative estimate for 

w := {ln$-u*) + (5.24) 

depending on the parameter p G [1, +oo). It is understood that the values of the constant c do 
not depend on p. We define 

u:=lni9, w n := min{n, (u — u*) + }, and v n := . (5.25) 

By arguing as done for the first estimate, we see that v n G L 2 (0, T; Vb), thus an admissible test 
function for (|2.28|) . and that Lemma 15. II can be applied with <3? = ip n given by (|5.14|) . by noting 
that v n = ip' n {§) and letting tp n (0) = 0. Therefore, as ?/> n (i?o) = by (I2.21|) . we obtain 

[M{W))+ I Vu-Vv n + [ W)-/3(P))v n = [ fv n (5.26) 
Jn JQt JQt JQt 

where / is still given by (|5,17p . Thanks to Lemma 15.41 we immediately derive that 

The next term on the left-hand side of (|5.26p is easily treated as follows 

/ Vu-Vv n = (2p-l) [ w 2 n V- 2 \Vw n \ 2 = [ |V<| 2 >-/ |V<| 2 

JQt JQt P JQt PJQt 

and the last one is nonnegative, since v n > and f3{'&) > /3(i?*) where v n > 0. In order to deal 
with the right-hand side, let us observe that / belongs to L q (Q) thanks to (|5.22[) ()5.23|) and 
(12, 13ft . Therefore, if q' denotes the conjugate exponent of q, we have 

fv n < \\f\\L«(Q)\\Vn\\ Lq ' iQ) < c||^ P_1 \\ Lq > (Q) = c|K||^ 2 ^ /p(Q) . 

Collecting the above estimates, we obtain 

MMWWh + f \V<?<cp\\<\\%£% for every te[0,T]. 
J Qt 

As both terms on the left-hand side are nonnegative, each of them satisfies the same bound. 
Therefore, owing to (|2.42|) . we derive that 

Il7/>P|| 2 < C filial! 2 _L ||V9/)P|| 2 o ^ < ™ H7/lP|l (2P ~ 1)/P 

W W n\\LW/i[Q) ^ c y\\ w n\\L°°(0,T;H) + II V W n \\L 2 {0,T;H) J - C P \\ w n ll ig '(2 P -i)/p(Q) • 

At this point, we note that (|5.24p and $ E L 2 (Q) trivially imply that w € L T (Q) for every 
r € [1, +oo). So we let n tend to infinity and conclude that 

ll?/)P|| 2 , < rnll7/; p ll (2p " 1)/p 

W W 1 1 1,10/3 (Q) ^ C P\\ W H L9 'C2 P -l)/p(Q) 

where w is given by (I5.24p . In other words, we have 

llwll m /, , < (m) 1 /^ lh/;ll (2p " 1)/(2p) 

IFIIil0f/3(Q) S (CP) \\ W \\ Lq '(2p-l)( Q y 

Finally, using the Holder inequality and terming \Q\ the Lebesgue measure of Q, we infer that 

ll-ll^/a^ < (c P )V(^| Q |V(^,o ii^ii^ 1 ^ 

< (c |Q|i/cw))V(*y/(*) M <wm < {cp) m P ) M <ww 
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since (Qj 1 /^ 9 ) is bounded with respect to p > 1. As we can assume the last constant c to 
be > 1, we conclude that 

\\ w \\L 10 P/ :i (Q) — \ C P) \\ w \\ L 2 Pq >(Q) \ b - 2 <) 

with c > 1. 

Conclusion of the proof. We rewrite (|5.27|) in the form 

IMIw W) < M^Mkllit;^ - : =^7 (5-28) 

and observe that a > 1 since </ < 5/3 by f)2.33l) . Now, we apply <|5.28|> to the divergent sequence 
{pk} defined by p k := a k and obtain 

\\ w \\ L *v k+iq > (Q) ^{cp k ) \\w\\ L2pkql(Q) . \p./M) 

Setting for convenience 

4 := m+ ll W ll L 2p fc g'(Q) 

and taking the positive part of the logarithm of both sides of (|5.29|) . we derive that 



4+1 < 7, — ln ( c Pfc) + h<7T~ ln ( c Pk) + 4 

%Pk 2p k 2p k 

the logarithms being nonnegative since c > 1. As this holds for every k > 0, we have that 

OO j 

ln + ||u;|| L2pfc9 / (Q) = 4 < C := 4 + ^ — ln(cpj), whence \\w\\ L 2 Pk«' (Q) - eC ( 5 - 30 ) 

i=0 ^' 

by noting that the series actually converges since pi = a 1 with a > 1. 

At this point, we can easily conclude the proof. Indeed, from (|5.30p . we immediately deduce 
that w G L°°(Q). Hence, coming back to (|5.24p . we derive that i? is bounded from above. 



6 Uniqueness 

In this section, we prove Theorem 12.71 The tool we use is the operator 31 : V * — > V$ given by 
the Riesz representation theorem, namely 



for v* G Vq and v G Vq, v = Jlv* means v* = -Av. (6.1) 



We note that 



-Ait, D£i>) = uv for every u G Vo and v £ H 
Jn 



(6.2) 



(u*,3to*) = (itV*)* for every it*, u* G F * (6.3) 



for every it G fl^O, T; V *) and for a.a. t G (0, T). (6.4) 

In (|6.4I) . || • ||* is the norm in V^* dual to the norm v *— » ||Vw||h in Vo, and the symbol on 
the right-hand side of (|6.3p is the corresponding inner product. By the Poincare inequality, 
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such norm and product in V * are equivalent to the standard ones and we mainly use them for 
convenience. Moreover, we recall that R = it — (3 satisfies assumption (|2.38p . Finally, despite of 
the general rule explained at the end of Section [21 we decide to compute all the constants we use 
in our estimates with care. In particular, we denote by Mq a constant that makes the following 
relations true 

||u||jf-i(n) — ^h|MI* for every v £ V * and \\v ||* < Mq||u||^ for every v G H (6.5) 

as well as the analogous of (|2.40|) for L 4 (f2). In (|6.5|) . the first norm is the standard one 
in V * = iJ-^O). 

To prove our uniqueness result, we have to show that any pair of solutions and 
(#2,^2, £2) to problem fl2~2"gj) - ([2~3"T|) satisfying the regularity requirements (^31 - 11X271) and 
having the 1? and X component bounded coincide, i.e., $1 = $2 and X\ = X2. So, pick such 
solutions. In particular, we can define M to be the maximum of the L°°-norms of the four 
functions #1, $2, and X2. We note that G, F', and G' are Lipschitz continuous on [— M, M], 
since they are smooth by (|2.3[) . In the sequel, L is the maximum of their Lipschitz constants on 
such an interval. Moreover, we can apply Remark 12.61 to (Xi and) X2 since $2 is bounded, and 
deduce that (1237]) holds for X 2 . 

Now, we write equations (|2.28p for both solutions and test their difference by — i^)- At 
the same time, we write equations (I2.29P for both solutions and test their difference by fidt(Xi — 
X2), where /1 £ (0,1) is a parameter whose value will be chosen later on. Finally, we sum 
the equalities we get to each other, rearrange, and add the same integral to both sides, for 
convenience. If we use the notation $ := $1 — $2 and X := X\ — X2, owing to (I6,2p - (|6.4p 
we obtain 

l\\m\\l+ [ CML-Ma)tf + / \dtX\ 2 + ^ [ |VX(t)| 2 + ^ / \x(t)\ 2 

z JQt JQt 1 JCl 1 JQ 

(d t G(X 1 (s))-d t G(X 2 (s)),ti(s))js+ [\(R(^)-R^2))(s),m)js 

J 

+ li I (F'(X 2 ) - F'(Xi)) t X + fx [ (G'(X2)$2-G\X 1 )$ l )d t X+^ I \X(t)\ 2 . (6.6) 
JQt JQt 1 Jn 

The only term on the left-hand side that needs some treatment is the second one. We have 

f (In #1 -1x102)0 > 4 f |0| 2 
JQt M JQt 

since < $j < M for i = 1,2. We deal with the first term on the right-hand side and use (|2.37|) 
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for X 2 as mentioned above. We get 

f\d t G{X 1 {s)) - d t G(X 2 (s)), ti(s))*ds 
Jo 

= [ {G'(X 1 (s))d t X(s) + (G'(X l (s))-G'(X 2 (s)))d t X 2 (s)^(s)) jf ds 
Jo 

<M Q [ WCiX^dtX^WnW^s^ds 
Jo 

+ M n [ ||G"(X 1 ( S ))-G / (X 2 ( S ))|| L 4 (n) ||,9 t X 2 ( S )|| i 4 (n) |K S )||^ S 
Jo 

< Mq L I \\d t X{s)\\ H \\&{s)\Uds + MnL ( \\X{s) || L 4 (n) ||$X 3 (s)|| L 4 (n) ||i^0ll.<k 
Jo Jo 

<£/ \d t x\ 2 + 2 -^_ fpwwld, 

° JQt V Jo 

ft u2r2 ft 

+ J o \\X(s)\\l m ds + SL- J \\d t X 2 (s)\\ 2 LHn) Ms)\\lds 
8 JQt V Jo 



f f Jyl-i J ^ ft 

+ Ml \VX\ 2 + M 2 |X| 2 + ^— / \\d t X 2 {s)\\l Hn) \\${s)\\lds. 
J Qt J Qt 4 Jo 

Moreover, we have 

[\(R(A)-R($2))(s),$(s))js<M n X R /W)klWs)IMs 
Jo Jo 

-1m I \^ 2 + mm ^r /VtollS* 

J J 

M / {F , (X 2 )-F , (X 1 ))a t X<^ f \d t X\ 2 + 2^L 2 I |X| 2 
JQt ° JQt JQt 

ii\ {G l (X 2 )$ 2 -G l (X 1 )$ 1 )d t X = fi [ {G\X 2 )-G'{X 1 ))^ 1 d t X- i x j G'(X 2 )tfd t X 



<iiLM [ \X\\dtX\+nL [ \0\\dtX\ 
JQt JQt 

/ \d t X\ 2 + 2^L 2 M 2 [ \X\ 2 + 2^L 2 [ 
4 J Qt JQt Jq, 

Finally, we treat the last term of (|6.6f) this way 



2 



/ \x(t)\ 2 = nf xd t x<^ I \d t x\ 2 + f , I |xp 

Jn JQ t ° J Qt JQt 



At this point, we observe that s 1— ► ||c^X 2 (s) || 2 4 ^ belongs to L 1 (0,T) by (|2.37j) . Hence, if 

we choose /1 in order that 2fiL 2 = 1/(4M), we see that (I6.6D and all the estimates we have 
performed allow us to apply the Gronwall lemma. This yields # = and X = a.e. in Q, whence 
the solutions coincide and the proof is complete. 
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